A Compact Difference Scheme for the Biharmonic Equation in Planar Irregular Domains

We present a finite difference scheme, applicable to general irregular planar domains, to approximate the biharmonic equation. The irregular domain is embedded in a Cartesian grid. In order to approximate ∆Φ at a grid point we interpolate the data on the (irregular) stencil by a polynomial of degree six. The finite difference scheme is ∆2QΦ(0, 0), where QΦ is the interpolation polynomial. The interpolation polynomial is not uniquely determined. We present a method to construct such an interpolation polynomial and prove that our construction is second order accurate. For a regular stencil, [7] shows that the proposed interpolation polynomial is fourth order accurate. We present some suitable numerical examples.